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ABSTRACT 

We have used the Submillimeter Array (SMA) to make the first interferometric observations (beam 
size ~ 1", or ~ 400 pc) of the 12 CO J=6-5 line and 435 ptm (690 GHz) continuum emission toward the 
central region (half power field of view 17") of the nearby ultra-luminous infrared galaxy (ULIRG) Arp 
220. These observations resolve the eastern and western nuclei from each other, in both the molecular 
line and dust continuum emission. At 435 fim, the peak intensity of the western nucleus is stronger 
than the eastern nucleus, and the difference in peak intensities is less than at longer wavelengths. 
Fitting a simple model to the dust emission observed between 1.3 mm and 435 fim suggests that 
dust emissivity power law index in the western nucleus is near unity and steeper in the eastern 
nucleus, about 2, and that the dust emission is optically thick at the shorter wavelength. Comparison 
with single dish measurements indicate that the interferometer observations are missing ~ 60% of 
the dust emission, most likely from a spatially extended component to which these observations are 
not sensitive. The 12 CO J=6-5 line observations clearly resolve kinematically the two nuclei. The 
distribution and kinematics of the 12 CO J=6-5 line appear to be very similar to lower J CO lies 
observed at similar resolution. Analysis of multiple 12 CO line intensities indicates that the molecular 
gas in both nuclei have similar excitation conditions, although the western nucleus is warmer and 
denser. The excitation conditions are similar to those found in other extreme environments, including 
the nearby starburst galaxy M82, the AGN hosting ULIRG Mrk 231, and the high-z QSO BR 1202- 
0725. Simultaneous lower resolution observations of the 12 CO, 13 CO, and C 18 J=2-l lines show that 
the 13 CO and C 18 lines have similar intensities, which suggests that both of these lines are optically 
thick, or possibly that extreme high mass star formation has produced in an overabundance of C 18 0. 
Subject headings: galaxies: individual (Arp 220), galaxies: ISM, galaxies: nuclei, galaxies: starburst 
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1. INTRODUCTION 

Galaxy-galaxy mergers are the phenomena violent 
enough to disturb the galactic potentials of both of the 
colliding galaxies. They create shocked gas at the merg- 
ing region, make interstellar matter fall into the new 
galactic potential, and induce strong starbursts. This 
extreme starburst heats up the surrounding dust, which 
sometimes radiates 10 12 L Q or more in infrared lumi- 
nosity. These sources are often called ultra-luminous in- 
frared galaxies (ULIRGs). 

Since molecular gas and dust create stars, observa- 
tions of these components toward ULIRGs are important 
to understand the nature of extreme starbursts. Mul- 
tiple molecular lines or multiple transition studies re- 
vealed that a large fraction of molecular gas in ULIRGs is 
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dominated by dense and warm gas (e.g., Solomo n et al.l 
119921 iGao fc Solomonl 12004 iPapadopoulos et all 120071 ) . 
In addition, the efficiency of star formation is tightly 
correlated with the dense gas fraction in m olecular gas 
dSolomon et al.lll99l IGao fc Solomo 3l200"l . 

ULIRGs are also important for the study of high-z sub- 
millimeter galaxies. The Submillimetre Common-User 
Bolometer Array (SCUBA) on the James Clerk Maxwell 
Telescope (JCMT) has d etected many high-z submillime- 
ter galaxies (SMGs; e.g .iHughes et al.lll998t[Eales et all 
[i"999t [20001 : iScott et al.l [20021 : iBorvs et al J 120031) . Most 
of them have infrared luminosities of > 10 12 L Q , and 
bright ones have ~ 10 13 Lq. Local ULIRGs can be stud- 
ied as nearby counterparts of high-z SMGs or used as 
nearby templates to derive photometric redshifts. To 
compare redshifted emission lines or continuum emission 
in high-z galaxies to those in local ULIRGs, we need 
to observe local ULIRGs at the same wavelengths as 
the rest wavelengths of the detected lines or continuum 
from the high-z galaxies. Many of the high-z galaxies 
were observed at mill i meter-wave (around 1-3 mm; 
e.g., iGreve et all 120051 : iTacconi et all 12006ft . Therefore 
the rest wavelengths of the detected lines or continuum 
are at submillimeter wavelengths. Our knowledge about 
submillimeter lines or continuum from local ULIRGs, es- 
pecially from compact starburst regions or from nuclei, 
is very limited so far, because only recently have high 
resolution submillimeter observations been possible. 

Here we present the first interferometric 12 CO (6-5) 
line and 435 /im (690 GHz) continuum observations of 
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Arp 2 20 using the Submillimeter Array (SMA; iHo et all 
[2001 . The interferometric 12 CO(2-l), 13 CO(2-l), and 
C 18 0(2-l) lines and 1.3 mm (226 GHz) continuum ob- 
servations have also been made simultaneously with the 
12 CO(6-5) and 435 /jm observations. Both 12 CO(6-5) 
and C 18 0(2-l) lines have never been observed (even in 
a single-dish telescope) toward this galaxy so far. Arp 
220 is the nearest ULIRG (79.9 Mpc, 1" = 387 pc; 
ISanders et al.ll2003h . and has therefore been well studied 
in various wavelengths. Optical images show distu rbed 
faint structures ("galaxies with adjacent loops"; lArd 
Il96fil) . which seem like re mnants of tidal tails pro duced 
by a galaxy-galaxy merger (j.Ioseph fc WrightJ[j~985h . sim- 
ilarto numerical simulation results (e.g.. lHernquistlll992L 
Il993fh The infrared lumino sity at 8 - 1000 p of 



this galaxy is 1.6 x 10 12 L Q (|Sanders et all 12003), an d 



high spatial r esolution radio contin uum (|NorrisT 19881 ). 
near-infrared (Graham et al.l 11990). and mid- infrared 
(jSoifer et al.l 119991 ) observations revealed two nuclei at 
the center with a separation of ~ 0.95" (<~ 370 pc). 
Molecular gas is very rich in this galaxy, ~ 9 x 10 9 M Q , 
and two-t hirds of this mass is concentrated within 400 pc 
in radius (jScoville et al.lfl997h . Such high gas mass con- 
centration is_amha£_alsj3_^ 

suits (e.g., iBarnes fc HernquistJll99ll : lMihos fc HernquistJ 
Il996h . These observations and numerical simulations in- 
dicate that Arp 220 is in the final stage of merging. 

The past high spatial resolution millimeter-wave inter- 
ferometric observations show molecular gas concentra- 
tions at the two nuclei, and these are embedded in an 
extended (~ 1 kpc) molecular structure, which seems to 
be a rot ating disk coinciden t with a dust lane in optica l 
images (|Scoville et al.lll997t IT5ownes fc Solomon! Il998f ). 
The molecular gas peak at each nucleus shows a steep 
velocity gradient with the direction of the gradient dif- 
ferent from that of the large-scale molecular gas disk. 
This suggests that a small-scale molecular gas disk ro- 
tates around each nucleus, and b oth disks are embedde d 
in the large-scale rotating disk (jSakamoto et al.l fl999f ). 
Recent ~ 0'.'3 molecular gas i mages resolved the de- 
tailed nuclear gas dist ributions (jDownes fe Eckaral2007t 
ISakamoto et alj 120081 ). which are consistent with the 
Hubble Space Telescope NICMOS near-infrared imaging 
results that suggest an op aque disk around one of the 
nuclei (jScoville et al.lll998l) . 

Dust emission als o peaks at the two nuclei. 
ISakamoto et alJ (|1999l ) reported that most of the contin- 
uum flux at 1.3 mm come s from the two nuclei, while 
iDownes fc Solomon! (|1998[ ) mentioned that half of the 
1.3 mm continuum flux comes from the extended com- 
ponent. It is suggested that a few tenths of the 8 60 
i|Sakamoto et al.ll2008D and 24.5 /im (jSoifer et al.lll999f) 
continuum flux comes from the extended component. 
Far-infrared observations with the Infrared Space Ob- 
servatory (ISO ) also suggest the existence of extended 
dust emission (|Gonzalez- Alfonso et a l. 2004). 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the center of Arp 220 with the SMA on 
March 2nd, 2005. The phase reference center was at 
a(2000) = 15 h 34 m 57!19 and 5(2000) = 23°30'11'.'3. The 
225 GHz atmospheric opacity was between 0.03 and 0.04, 
which was measured at the nearby Caltech Submillimeter 
Observatory. Six of the eight 6 m antennas were used 



with projected antenna separations between 14 m and 
68 m. High frequency receivers were tuned to observe the 
redshifted 12 CO(6-5) line (679.13 GHz) in the lower side 
band (LSB) and the 689.13 GHz (~435 fim) continuum 
emission in the upper side band (USB). Low frequency 
receivers were tuned to observe the redshifted 13 CO(2-l) 
line (216.47 GHz) and C 18 0(2-l) line (215.64 GHz) in 
the LSB and the redshifted 12 CO(2-l) line (226.42 GHz) 
in the USB. The double sideband system temperature for 
the high frequency band ranged from 2000 K to 2500 K 
for most of the time (i.e., at high elevation), and that for 
the low frequency band from 140 K to 180 K. The SMA 
correlator covers a 2 GHz bandwidth for each sideband of 
both high and low frequency bands, which corresponds 
to the velocity range of ~ 880 km s _1 and ~ 2700 km s _1 
for each sideband of the high and low frequency bands, 
respectively. The channel width was configure to have 
3.25 MHz (~ 1.4 km s^ 1 ) and 0.8125 MHz (~ 1.1 km s" 1 ) 
for the high and low frequency bands, respectively. 

We calibrated the data using the Owens Valley Ra- 
dio Observatory software package MIR, which is modi- 
fied for the SMA. For the high frequency band calibra- 
tion, we used a partially resolved source, Callisto, as a 
gain (amplitude and phase) and flux calibrator 9 , since 
quasars were too weak for gain calibration. Callisto was 
about 50° away from Arp 220. The r.m.s. phase fluc- 
tuation was about 21°. We used a source model at the 
gain calibration to correct for the effect of the partially 
resolved structure. Bandpass calibration was done us- 
ing three sources, Mars, Callisto, and Ganymede to gain 
the signal-to- noise ratio (S/N). Ceres was imaged after 
the calibrations, and its flux was 22% lower than the 
calculated flux 9 . Although the flux error for Ceres was 
about 20%, we conservatively adopt the flux error of 30% 
for the high frequency band data throughout this paper. 
For the low frequency band calibration, we used Ceres, 
which was about 40° away, as the gain calibrator due 
to its closeness, and Mars, Callisto, and Ganymede were 
used as bandpass calibrators. Callisto was imaged after 
the calibrations, and showed 21% lower than the calcu- 
lated flux. Hereafter, we adopt the flux error of 20% for 
the low frequency band data. 

Data from 5 antennas were used for the 12 CO(6-5) line 
imaging because of a correlator problem. For 435 //m 
continuum imaging, we used the data from all 6 anten- 
nas after discarding the problematic data. We subtracted 
the continuum emission from the line emission data be- 
fore line imaging. Since the 12 CO (6-5) line width is com- 
parable with the band width, we could not obtain the 
continuum emission from the same sideband. We there- 
fore subtracted the continuum emission using the other 
sideband (i.e., USB) in the uv plane using the National 
Radio Astronomy Observatory software package AIPS. 
For the low frequency (1 mm) data, we created the con- 
tinuum image using the line free channels in the same 
band. The line images for the low frequency data were 
made after the subtraction of the continuum emission 
from the data in the uv plane. 

The calibrated data were binned, and the final chan- 
nel maps have velocity resolutions of about 30 km s _1 

9 The flux values for Callisto and Ceres are inferred 
from the SMA Planetary Visibility Function Calculator 
(http: //smal . sma.hawaii . edu/planetvis . html ) 



SMA 12 CO(J = 6-5) and 435 /im observations of Arp 220 



TABLE 1 

Parameters for the continuum and molecular line images 



Wavelength (Frequency) 


Synthesized Beam Size 


Velocity Resolution 


R.M.S. Noise 


or Line 


and Position Angle 






[fim (GHz)] 


(Linear scale) 


[km s~ ^] 


[mjy beam" - ""'" (niK)] 


435 (689.13) 


T/2 x 0'.'9, 139° 


— 


190 (450) 




(470 pc x 350 pc) 






1320 (226.46) 


3'.'8 x 3'.'3, 28° 




4.8 (9.1) 




(1.47 kpc x 1.28 kpc) 






1380 (216.46) 


3'.'9 x 3" 5, 27° 




5.7 (11) 




(1.51 kpc x 1.36 kpc) 






"CO(6-5) 


T/3 x O'.'S, 129° 


30.1 


535 (1400) 


12 CO(2-l) 


(500 pc x 310 pc) 






3('8 x 3'.' 3, 28° 


5.4 


32.2 (61) 


(1.47 kpc x 1.28 kpc) 




13 CO(2-l), C ls O(2-l) 


3'.'9 x 3" 5, 27° 


30.8 


13.1 (25) 




(1.51 kpc x 1.36 kpc) 







for the 12 CO(6-5), 13 CO(2-l), and C 18 0(2-l) lines, and 
about 5 km s -1 for the 12 CO(2-l) line. We use the radio 
definition for the LSR velocity in this paper, which is 
u lsr = c(l — ^/t'rost)- The C 18 0(2-l) line was observed 
up to the LSR velocity of 5550 km s -1 , covering about 
72% of the total line width, since this line was at an edge 
of the bandpass. 

We CLEANed the images with natural weighting, and 
the resulting synthesized beam sizes were about 1" for 
690 GHz band images, and about 3" - 4" for 230 GHz 
band images. The beam sizes and the r.m.s. noise lev- 
els for all the images are summarized in Table [TJ The 
half-power width of the primary beam at 690 GHz and 
230 GHz are 17" (6.6 kpc) and 52" (20.1 kpc), respec- 
tively. These sizes are much larger than the sizes of the 
line/continuum emitting regions (at most a few arcsec- 
onds), and we did not make any primary beam correction 
to our images. 

Our 12 CO (6-5) images show a systematic position shift 
of about 0"7 from the peak positions of the double nu- 
cleus reported in the previous observ ations at longer 
wavelengths (e.g.. lSakamoto et al.lll999h . This offset can 
be explained by the baseline error, which is about 0.3A 
or less at 690 GHz. Furthermore, we used the phase cal- 
ibrator with about 50° away from the source. Hence we 
shifted the positions of the 12 CO(6-5) line images based 
on the peak positions and kinematics information to be 
consistent with the previously published results. The 
435 fim continuum image also shows about Qf.'l position 
shift, which is again explained by the baseline error, and 
therefore shifted the position with the same manner as 
the 12 CO(6-5) data. The low frequency band images did 
not show any noticeable position shifts, and therefore we 
did not shift any images. 

3. RESULTS 

3.1. 435 \xm Continuum Emission 

The 435 /im (689 GHz) continuum emission image is 
shown in Fig.[T] The image clearly shows two peaks with 
a separation of about 1", consistent with the past con- 
tinuum images in centimeter, millimeter, and infrared 
wavelengths. We therefore call these peaks eastern and 
western nuclei as is the past studies. Intensities of the 
western and eastern nuclei are 1.28 ± 0.38 Jy beam -1 
and 0.96 ± 0.29 Jy beam -1 , respectively. The total flux 
density of the continuum emission is 2.5 ± 0.8 Jy. Since 
the flux distribution can be smeared by the phase fluctu- 
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Fig. 1. — The 435 (im (689 GHz) continuum emission im- 
age of the central region of Arp 220. The synthesized beam 
(T/2 x 0'.'9 or 470 pc X 350 pc) with the P.A. of 139° is shown at the 
botto m-left corner. The two crosses indicate the 1.3 mm continuum 
peaks (Sakamoto ct al. 1999; Downcs & Solomon 1998), which cor- 
responds to the double nucleus. The contour levels are 3, 4, 5, and 
6a, where la = 190 mjy beam -1 . The position offsets are mea- 
sured from a(2000) = 15 h 34 m 57!25 and (5(2000) = 23°30'lT/4. 



ation or baseline errors, we convolved the image to larger 
beam sizes (up to 10" with a Gaussian convolution) and 
measured the total flux density. It did not change from 
2.5 Jy, indicating that the flux smearing effect due to the 
phase or baseline errors is small. The r.m.s. phase fluc- 
tuation at this wavelength was indeed only 21° (Sect. [J), 
which induces < 0"1 smearing, so it is consistent with 
this result. 

The previously published 450 fj,m (670 GHz) sin- 
gle dish results show a very large variation in the de- 
tected flux densities; the United Kingdom Infrared Tele- 
scope (UKLR T) UKT14 result shows the flux density of 
3.0 ± 1.1 Jy (|Eales et al.lll989t ). but th at of the JCMT 
SCUB A result shows 6.286 ± 0.786 Jy (|Dunne fc Ealed 
l2001h . We plot the submillimeter (1.5 mm - 300 /im 
or 200 GHz - 1000 GHz) spectrum energy distribution 
(SED) using published single dish results in Fig. [2] The 
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200 400 600 800 1000 

Frequency [GHz] 

Fig. 2. — Submillimeter wavelength (1.5 mm — 300 fim 
or 200 GHz - 1000 GHz) SED for various components in 
Arp 220. Crosses, squares, diamonds, and pluses indicate the 
continuum flux density of the entire system, the eastern nu- 
cleus, the western nucleus, and the extended component, re- 
spectively. A down-arrow indicates the 3a upper limit for the 
extended component. Bars on the observational points indi- 
cate the observational errors. Solid, dotted, dashed, and dot- 
dashed lines are the \ 2 fitting for each component. The low fre- 
quency data for the double nucleus a nd th e exte nded component 
are ta ken from IDownes fe Solo mon! (1998) and Sakamot o et al.l 
(fT999T) for 229 GHz IWiedner et al.l | |2002T) for 343 GHz, and 
ISakamoto et al.l 1 120081) for 345 GHz. The data for the total flux 
densit y are taken fro m following pa p ers: 2 14 GHz: Woo dv~et al.l 
1 19891) . 229 GHz: IScoville et al.l (fl99W l; IDownes fc Sqlomqnl 
| 1S?98I) : ISatemoto et al.l (I1999I). "240 GHz: ICarico et al.l (119921), 
273 G Hz: IRigopoulou et al.l {19961). 3 53 GHz: IDunne fc Ealesl 
lUSHD, 3 75 GHz: |Eales et al.l (119891): IRigopoulou et al.l (119961 )' 
667 GHz: lEales et al.l U989l);IDunne"& ; Eales ( 200T|), and 857 GHz: 
lEales et al.l (119891) : IRigopoulou et al.l 1 1199611 

large crosses in the plot are the data for the total flux 
densities at various frequencies, and the solid line shows 
the x 2 fitting of the data (see Sect. 14. ll for details). As 
shown in the figure, the JCMT SCUBA data point is 
on the fit, but the UKIRT UKT14 data point is signifi- 
cantly lower than the fit. There is no report for the time 
variation of flux density at submillimeter wavelength in 
this galaxy so far. In addition, the measurements for the 
total flux density in the plot is distributed over three 
decades, but only 450 /im shows the large difference, 
which is unlikely to be caused by the time variation. We 
therefore assume that the JCMT SCUBA value is more 
accurate. Using the fit, we estimated the total flux den- 
sity at 435 fim as 5.9 Jy. This suggests that our 435 ^m 
continuum observations missed ~ 58% of the total flux 
density. 

Since this missing flux is larger than our flux error of 
~ 30%, and the effect of the flux smearing due to phase 
or baseline errors is small, this missing flux is probably 
due to the existence of an extended component. If the 
missing flux of 3.4 Jy is due to an extended component 
with a Gaussian distribution that has a full width at 
half maximum of ~ 3", its peak will be only twice of 
our r.m.s. noise. We did not detect any significant sig- 
nal with larger beam as mentioned above, but the r.m.s. 
noise level increased in the larger beam images due to a 
small number of data points at shorter baselines. It is 
therefore not clear whether the missing flux in our data 
is due to the lack of shortest baselines or to low S/N. 

3.2. 12 CO(6-5) Line Emission 

The channel maps of the 12 CO(6-5) line emission at the 
nuclear region of Arp 220 are shown in Fig. [3l and the 
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Fig. 3. — Channel maps of the 12 CO(6-5) line emission. Con- 
tinuum emission is already subtracted. The LSR velocity (radio 
definition) in km s" 1 is shown at the upper-left corner of each 
channel map, and the synthesized beam (1"3 X 0''8 or 500 pc 
X 310 pc) with the P.A. of 129° is shown at the lower-left cor- 
ner of the first channel map. The position offsets are measured 
from a(2000) = 15 h 34 m 57?25 and 5(2000) = 23 o 30'lr74. The two 
crosses in each channe l map indicate the 1.3 mm continuum peaks 
(Sakamoto et al. 1999). The contour levels arc —3,3,4,5, and 6<r, 
where la = 535 mjy beam -1 (= 1.4 K). 

integrated intensity (moment 0) and intensity-weighted 
mean velocity (moment 1) maps are shown in Fig.[3J The 
channel maps reveal that the molecular gas traced by 
the 12 CO(6-5) line exhibits different kinematics around 
each nucleus; molecular gas associated with the east- 
ern nucleus tends to display gas kinematics moving from 
south-west to north-east, and that associated with the 
western nucleus from south-east to north-west. These 
features can also be seen in the moment 1 map for the 
eastern nucleus. Molecular gas around the western nu- 
cleus in the moment 1 map, on the other hand, does 
not exhibit a clear velocity gradient as in the channel 
maps. This may be due to the nature of the moment 
1 map (i.e., intensity weighted velocity field map) with 
weak emission at the northwestern region. The kine- 
matics on spatial scales larger than that immediately 
associated with the individual nucleus shows a velocity 
gradient along north-east to south-west direction with 
much shallower gradient than that of the eastern nucleus. 
All of these kinematics features mentioned above are 
also seen in the p revio usly pu blished lower- J 12 CO lines 
(ISakamoto et alJ fl99flL 120081 E owncs fc Solomon! I~1998: 
IScoville et al.lll997fl . The kinematics in the west nucleus 
displays a smaller velocity gradient in our data, prob- 
ably due to poorer spatial resolution and lower S/N of 
our data. The detection of the large-scale kinematics 
feature suggests that our data detected some of the ex- 
tended (a few arcsecond scale) component, but much less 
than that detect ed in 12 CO(1-0), 12 CO(2-l), or 12 CO(3- 
2) observations (IScoville et al.lll997t IDownes fc Solomon! 
1998; ISakamoto et al.lll999L l2008fK 
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Fig. 4. — (a) Integrated intensity (moment 0) and (b) intensity- 
weighted mean velocity (moment 1) maps of the 12 CO(6-5) line 
emission. Continuum emission is already subtracted. The po- 
sition offsets, the synthesized beam, and the crosses are the 
same as in Fig. \3\ The contour levels for the moment map 
are (2, 4, 6, 8, 10, 12 and 14) X 68 Jy beam -1 km s -1 (= 
175 K km s — 1 ). The contour levels for the moment 1 map 
are -90,-60,-30,0,30,60,90,120,150, and 180 km s" 1 , where 
km s — 1 corresponds to the LSR velocity of 5351 km s — 1 . Zero 
velocity is shown in thick solid contour, and negative and positive 
velocities are shown in dashed and solid contours, respectively. 

The integrated intensity image, on the other hand, dis- 
plays a single peak that is elongated along the two nuclei, 
which looks different from our 435 /im continuum im- 
age or the previously publi shed high resolution 12 CO (2- 
1) and 12 CO(3-2) images (|Sakamoto et al.l I1999L 120081 : 
iDownes fc Solomonlll998l ) at first glance. Our image is 
rath er similar to the low resolution maps of 1 2 C O ( 2- 1 ) 
line (jSakamoto et alJll999t IScoville et al.lll997l ) without 
diffuse and extended (larger than a few arcsecond) emis- 
sion from their maps. The total integrated 12 CO(6-5) in- 
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Fig. 5. — Spectra of the 12 CO(6-5) line at the two nuclei. 
Continuum emission is already subtracted. 

tensity is 1250 ± 250 Jy km s _1 in our data. Since there 
is no published result for the 12 CO(6-5) line emission, 
we could not estimate the missing flux of our 12 CO (6-5) 
data. 

Fig. [5] shows the 12 CO(6-5) line spectrum at each nu- 
cleus. The peak brightness temperature of the western 
nucleus is 8.7 ± 2.6 K (3.4 ± 1.0 Jy beam -1 ) around the 
LSR velocity of ~ 5300 km s _1 , and that of the east- 
ern nucleus 6.1 ± 1.8 K (2.4 ± 0.7 Jy beam -1 ) around 
~ 5500 km s -1 . The overall velocity ranges for the 
two nuclei are similar to those in the high resolution 
12 CO (2-1) and 12 CO(3-2) observations (jSakamoto et all 

3.3. 1.32 mm and 1.38 mm Continuum Emission 

Both 1.32 mm (226.46 GHz) and 1.38 mm 
(216.46 GHz) continuum emissions are detected at 
a significant signal level (~ 30er), and the emission is 
unresolved at both wavelengths at our resolution of 
3" — 4" (the 1.32 mm image is shown in Fig. [5^,; the 
1.38 mm image is almost the same as the 1.32 mm 
image, and not shown here). The total flux densities 
for 1.32 mm and 1.38 mm are 167 ± 33 mJy and 
160 ± 32 mJy, respectively. From the dust continuum 
fitting (Sect. 13. ip . we estimate the total flux densities 
at 1.32 mm and 1.38 mm as 153 mJy and 178 mJy 
(including the non-thermal flux contribution). The 
observed total flux density therefore agrees with that 
from the SED fitting within our calibration error. We 
therefore conclude that our 1.32 mm and 1.38 mm 
continuum data have no missing flux. 

3.4. 12 CO, 13 CO, and C 18 J=2-l Line Emissions 

The 12 CO(2-l) line image exhibits an extended struc- 
ture along north-east and south-west direction even with 
our low resolution image (Fig. [5b), which is consis- 
tent with the past interferometric maps (ISakamoto et aTl 
ll999HDownes fc Solomonlll998l : IScoville et al.lll997f) . On 
the other hand, the 13 CO(2-l) and C 18 0(2-l) line images 
(Fig. [5J:, d) are unresolved at our resolution. 

The integrated intensities of the 12 CO, 13 CO, and 
C 18 J=2-l lines are 1430 ± 290 Jy km s -1 , 45.7 ± 
9.1 Jy km s , and 31.5 ± 6.3 Jy km s , respectively. 
Our observation covers only ~ 72% of the line width of 
the C 18 0(2-l) line (Sect. [J), so the derived value should 
be a lower limit (see also Sect. I3.5|) . 

We compared the integrated intensities with the single 
dish results for the 12 CO(2-l) and 13 CO(2-l) lines. The 
JCMT observations of the 12 CO(2-l) line (21" -22" res- 
olutio n) indicate 1730 ± 350 Jy km s -1 (IWiedner et al.l 
l200l and 1549±311 Jykms -1 (|Greve et al.ll2008h . The 
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Fig. 6. — Continuum and line images taken in the low frequency 
band. Continuum is already subtracted from the line images. Syn- 
thesized beams are shown at the lower-left corner of each image, 
and their sizes are in Table[T] The crosses are the same as in Fig. [3] 
The position offsets are measured from q(2000) = 15 h 34 m 57!24 
and 5(2000) = 23°30'11'.'3. (a) 1.32 mm continuum. Contour lev- 
els are -3, 3, 5, 10, 20, and 30cr, where la = 4.8 mjy. (b) 12 CO(2-l) 
integrated intensity. The contour levels are 5, 10, 20, 50, 100, 150, 
350 X 2.5 Jy beam" 1 km s" 1 (= 4.8 K km s" 1 ). (c) 13 CO(2-l) 
integrated intensity. The contour levels are 3, 5, 10, 15, 20, and 
25 X 2.1 Jy beam" 1 km s" 1 (= 4.0 K km s" 1 ). (d) C ls O(2-l) 
integrated intensity. The contour levels are the same as in the 
13 CO(2-l) map. 

flux differences between our and these observations are 
within the calibration errors. The 13 CO(2-l) line was 
observed with the Institut de Radio Astronomie Mil- 
limetrique (IRAM) 30 m telescope (11" r esolution) and 
JCMT (21" resolution) (jGreve et al.ll2008l ). and they ob- 
tained the integrated intensities of 60 ± 13 Jy km s _1 
and 70 ± 16 Jy km s _1 , respectively. The differences be- 
tween our and these observations are again explained by 
the calibration errors. We therefore conclude that our 
12 CO(2-l) and 13 CO(2-l) line data have no missing flux. 

Since there is no single dish C 18 0(2-l) line observation, 
we cannot estimate the missing flux for this line. On the 
other hand, since we did not see any significant missing 
flux in the 1.32 mm and 1.38 mm continuum data and 
the 12 CO(2-l) and 13 CO(2-l) line data, we expect no 
significant missing flux in the C 18 0(2-l) line data. 

3.5. Line Spectra and Ratios 

To compare the line spectra and intensities of all the 
four obtained lines, we convolved the data into the largest 
synthesized beam size of 3" 9 x 3" 5 with the RA. of 27°, 
which is the beam size of the 13 CO(2-l) and C 18 0(2-l) 
lines. Fig. [7] shows the spectra of all these four lines. 
Note that the sensitivity toward extended structure (i.e., 
uv coverage) is different between the J = 2 - 1 lines and 
the J = 6 - 5 line in this figure. 

The 12 CO (2-1) line shows a double peak profile with 
stronger intensity at lower velocity, whic h is consistent 
with t he past single dish m easurements ( Wicdnc r et al.l 
I2002t ISolomon et ail Il990f ). The line width of the 
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Fig. 7. — Multiple CO line spectra toward the central region 
of Arp 220 at 3'.'9 X 3'.'5 resolution (P.A. = 27°). Thick solid, 
thin solid, thin dashed, and thin dash-dot lines are the spectra 
of the 12 CO(2-l), 13 CO(2-l), C ls O(2-l), and 12 CO(6-5) lines, re- 
spectively. Vertical axis is brightness temperature in Kelvin, and 
horizontal axis is LSR velocity in km s _1 . Due to the weakness of 
the 13 CO(2-l) and C ls O(2-l) lines, these intensities are increased 
by a factor of 10 in this figure. Since the C ls O(2-l) line is located 
at the edge of the bandpass, the spectrum finishes around the LSR 
velocity of 5550 km s _1 . Note that the uv coverage is different 
between J = 2 - 1 lines and the J = 6 - 5 line in this figure. 

12 CO(2-l) line reaches 900 km s" 1 (4900 - 5800 km s" 1 ). 
The 13 CO(2-l) and C 18 0(2-l) lines exhibit very similar 
line profiles and intensities. The 12 CO(6-5) line is mostly 
emitted at higher velocities and is weak at lower veloci- 
ties. This asymmetry is s imilar to that of the HCN(4-3) 
line (jWiedner et al.ll2002fl . 

We matched the shortest uv distance between the 
12 CO(6-5) and 12 CO(2-l) data to match the spatial 
structures between these two lines, and measured the 
12 CO(6-5)/(2-l) intensity ratio to be 0.34 ± 0.12. This 
value is much lower than the lower- J ratios, such as 
12 CO(3-2)/(2-l) of 0.85 ± 0.2 4 dWiedner et al. 1120021) or 
12 CO(2-1)/(1-0) of 0.65 ±0.1 (|Scoville et al.lll997lh This 
is mostly due to the different line profile of the 12 CO (6- 
5) line from the lower-J CO lines. This indicates that 
the lower velocity is dominated by lower-J CO lines, but 
the higher velocity is rich in highe r- J lines. This t r end is 
consistent with the observations of lWiedner et al.l (|2002D 
that the 12 CO (3-2) line intensity decreased relative to the 
12 CO(2-l) line in the lower velocity part, but stay almost 
constant at higher velocities. The relation between the 
spatial distribution and the velocities of the molecular 
gas in this galaxy is, however, not simple; all the molec- 
ular gas components in this galaxy, namely the two nuclei 
and the extended co mponent, have low and high veloc- 
ities (Fig. [5j see also ISakamoto et al.lfl999l l2008h . It is 
therefore difficult to tell which component contributes to 
high or low excitation condition only from the large-scale 
line spectra. The 12 CO transition ratios for each nucleus 
is derived in Sect. 14.31 

The 12 CO(2-l)/ 13 CO(2-l) line ratio is 13.0±3.7 at our 
resolution of 3"9 x 3'.'5. We also convolved our data to 
13" resolution (« single dish resolution), and the ratio 
was 16.2 ± 4.6. These values are similar or slightly lower 
than the values observed in U/LIRGs, and similar or 
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Fig. 8. — Continuum flux ratio between the eastern and western 
nuclei as a function of frequency. Our data point is the right-most 
one, and the data for 229 GHz and 343 GHz are the same as in 
Fig.0 

somewhat higher than starburst or Seyfert galaxies; the 
12 CO(2-l)/ 13 CO(2-l) line ratios in U/LIRGs observed 
with single dish telescopes ar e 23.5±4 (|Casoli et al.lfl 992') 
and 16±8 (Gl enn &: Hunterll2001[ ) , and those in star burst 
and Se yfert galaxies are 13 ± 5 (Aalto et al.l fl995f) and 
13 ± 1 (|Papadopoulos fc Seaauist 119981 ). respectively. 

The line profiles and the line intensities of the 13 CO(2- 
1) and C 18 0(2-l) are almost identical (Fig. O, and the 
13 CO(2-l)/C 18 0(2-l) line intensity ratio is 1.0 ± 0.3 be- 
tween the velocity range of 4800 km s _1 and 5550 km s _1 , 
which is the velocity range we observed the C 18 0(2-l) 
line (Sect. E2). The 13 CO(2-l) line has 92% of its to- 
tal integrated intensity in this velocity range, so the 
result with the full line width will not change signif- 
icantly. This very low 13 CO(2-l)/C 18 0(2-l) intensity 
ratio mat ches very well wit h 13 CO(l-0)/C 18 O(l-0) = 
1.0 ± 0.3 (jGreve et al.ll2008f l We discuss this ratio in 
Sect, tfjl 



4. DISCUSSION 

4.1. Dust Emissivity Index and Dust Opacity of the 
Two Nuclei 

We plot in Fig. [5] the flux ratios between the two nu- 
clei as a function of frequency. Although our continuum 
data have ~ 30% of calibration errors and a large amount 
(~ 60%) of missing flux, the flux density ratio between 
the two nuclei is accurate, since it depends on the noise 
level of the map. Hence this flux ratio diagram has higher 
accuracy than the SED diagrams shown in Fig. [2l The 
flux ratio in our data is 0.75±0.19. If we assume that the 
missing flux of our data is due to the extended emission 
(see Sect. I3.1| l . the size of the extended emission is larger 
than the separation of the two nuclei. Therefore the cor- 
rection of the missing flux will increase both the fluxes 
of the two nuclei with almost the same amount, and the 
flux ratio increases toward unity namely the ratio de- 
rived above is the lower limit. Contamination of the flux 
from one nucleus to the other due to the large beam size 
is small; this effect lowers the ratio by ~ 6%, which is 
smaller than the calibration errors. As frequency goes 
down to ~ 200 GHz, the flux ratio also goes down, and 
the flux of the eastern nucleus is about half or less of that 
of the western nucleus. Since the data for the two lower 
frequencies may not have s ignificant missing flux (e.g., 
I Sakamoto et al.l I1999L [2008) , the difference between the 



ratio at 689 GHz and those at lower frequencies will be 
more pronounced if we correct for the missing flux of the 
689 GHz data. Since the emission at these frequencies 
is dominated by dust emission, this result suggests that 
the dust SEDs are different between the two nuclei. 

We then made continuum SEDs for the two nuclei at 
submillimeter wavelength (Fig. \2§ as well as that for the 
total flux density to discuss the dust property differences 
more quantitatively. The crosses, filled squares, and filled 
diamonds in the plot are the data for the total flux den- 
sity, east nucleus, and west nucleus, respectively. The 
solid, dashed, and dotted lines are the \ 2 fitting of the 
data with a function eB(Td), where B(Td) is the Planck 
function for temperature Td and e is the emissivity func- 
tion. We adopted a form e — 1 — exp[— {y/v^\; v c is 
the critical frequency where the opacity is unity, and (3 
is the emissivity index. We adopted the source sizes of 
0'.'27 x 0'.'14 for the e astern nucleus and 0"1 6 x 0"13 for 
the western nucleus (|Sakamoto et al.ll"2008f ). since their 
images have the highest spatial resolution at the high- 
est frequency around submillimeter wavelengths. For 
the source size of the total flux density, we adopted the 
deconvolve d size of the extende d 12 CO(2-l) emission of 
l'/94x l'/28 (|Scoville et al.lH99l . The flux density for the 
non-thermal component is subtracted from all the fluxes 
in the plo ts before the f i tting, following the method ex- 
plained in lScoville et al.l (|1991h . although it is not signifi- 
cant especially at high frequencies; about several milli-Jy 
or several percent of the flux around 200 — 300 GHz, and 
even smaller at 600 — 700 GHz. We also estimated the 
CO flux contamination into the total flux, but most of 
the data are affected for only a few %, and the fitting re- 
sult did not change. The fitting results are summarized 
in Tabled 

We obtained from the fit to the total flux density Td of 
51 K and (3 of 1.4, which are consistent wit h the past esti- 
mate s ofT ri ~ 40-60K with/3 ~ 1.2-2.0 (IScoville et all 
19911: ID owncs fc Solomon! fl998t iDunne fc EaleT ^OOl; 



Klass et al.l |2001| )" The critical frequency is estimated 
to be 2200 GHz (~ 140 /mi), which is also consistent 
with the esti mation that dust is already optically t hick 
at - 100 ^m (IScoville et al.lH99ll iKlass et al.ll200l . If 
we adopt a smaller source size for the total flux of l'/13 
(|Scoville et al.lll997h . T d , 0, and v c would be 66 K, 1.3, 
and 2000 GHz (~ 150 /im), respectively. The tempera- 
ture is slightly higher, but the critical frequency and the 
emissivity index are still consistent with the past esti- 
mates. 

The fitting for the two nuclei gives high temperatures 
of 83 K and 180 K for the east and west nuclei, respec- 
tively. These high temperatures are due to the small 
source sizes. The emissivity indices are different between 
the two nuclei, 2.1 for the eastern nucleus and 1.2 for the 
western nucleus. The critical frequencies are estimated 
to be 370 GHz (~ 810 /im) for the eastern nucleus and 
190 GHz (~ 1.6 mm) for the western nucleus. 

The large uncertainties in our fitting for the two nuclei 
are the source flux density and the adopted source size. 
If we increase the 689 GHz flux density of each nucleus 
by 30%, which corresponds to the calibration error of our 
data may have (Sect.[2|), the fitting results would be Td = 
120 K and 310 K, (3 = 1.8 and 0.7, and v c = 520 GHz 
(~ 580 /im) and 610 GHz (~ 490 /im) for the eastern 
and western nuclei, respectively. We adopted the source 



Matsushita et al. 



TABLE 2 

Fitting results of dust spectral energy distributions 



Source 


Dust Temperature T<j 
[K] 


Emissivity Index (3 


Critical Frequency (Wavelength) v c 
[GHz (jjem)] 


Arp 220 (total) 


51 - 66 


1.3 - 1.4 


2000 - 2200 (~ 140 - 150) 


East Nucleus 


49 - 120 


1.8 - 2.1 


370 - 520 (~ 580 - 810) 


West Nucleus 


97 - 310 


0.7 - 1.2 


190 - 610 (~ 490 - 1600) 


Extended Component 


~ 38 


~ 2.4 


~ 1500 (~ 200) 



size derived at 860 fim ()Sakamoto et al.l [2008T ). but the 
effective source size at 435 /im may be different from that 
at 860 /tin. This is because the opacity is wavelength 
dependent, and therefore the effective source size is also 
wavelength dependent. Ideally we need to determine the 
source size at each wavelength, but our spatial resolution 
at 435 /xm is too low for this. Since opacity is higher at 
shorter wavelength, the effective source size at 435 fim. 
may be larger than what we adopted. If we increase the 
source size (area) of each nucleus by 50%, which roughly 
corresponds to the deconvolved size of the double nucleus 
in our data, the fitting results would be Td = 49 K and 
97 K, (3 = 2.1 and 1.1, and v c = 400 GHz (~ 750/xm) 
and 210 GHz (~ 1.4 mm) for the eastern and western 
nuclei, respectively. 

The fitting results, with the uncertainties, indicate that 
the dust temperature, the emissivity index, and the crit- 
ical frequency for the eastern nucleus are better con- 
strained than the past observations, because of our high 
frequency observations with high spatial resolution. The 
eastern nucleus seems to have a warm temperature of 
49 - 120 K, a steep emissivity index of ~ 2, and be- 
comes optically thick at frequencies above ~ 400 GHz. 
iSakamoto et alj (|2008l ) derived a dust temperature of 30 
- 160 K, so our estimate narrows the range. They also 
estimated the 350 GHz opacity of 2.8 for (3 = 2 (and 0.8 
for (3 = 1), hence our result is somewhat lower, but still 
both results indicate a high opacity condition at submil- 
limeter wavelengths in the eastern nucleus. 

The emissivity index and the critical frequency for the 
western nucleus are less constrained than those for the 
eastern nucleus, but are better constrained than the past 
observations. Our fitting results indicate a shallow emis- 
sivity index of about unity and a l ow cr itical frequency 
of < 600 GHz. iDownes fc EckartJ (120071) estima ted the 
230 GHz opacity as > 0.7 and ISakamoto et ail (|2008[ ) 
estimated the 350 GHz opacity as 0.8 — 5.3 for (3 — 2 
(the estimated ranges depends on the source size and 
the flux errors). Both of these observations and our re- 
sults indicate the western nucleus is optically thick at 
submillimeter wavelengths. 

The fitted temperature for the western nucleus, on the 
other hand, has a large range, and therefore it does not 
set a tighter limit than the past estimations. It appears 
that the western nucleus is warmer than the eastern nu- 
cleus, with dust t e mpera tures of a few hundred Kelvin. 
IDownes fc Eckartl (|2007l) suggested a dust temperature 
of 170 K, and ISakamoto et all (|2008l ) derived a temper- 
ature of 90 - 180 K, hence our estimate is consistent. 

Our results indicate that the derived properties, es- 
pecially the emissivity, indices are different between the 
two nuclei, suggesting that the dust properties, such as 
dust size distributions or dust compositions, are differ- 
ent. This difference may reflect the dust properties of the 



original host galaxies of each nucleus, or the difference in 
activities, such as star formation or AGN activities, in 
the nuclei after the merger. 

4.2. Extended Component in the Dust Emission 

Our 435 /xm (689 GHz) continuum data missed a sig- 
nificant amount (~ 61%) of the total flux. Here we dis- 
cuss whether this missing flux can be attributed to an 
extended component in the dust emission. The molec- 
ular gas clearly has an extended component with a size 
of ~ 2" (~ 1 kpc), which is inter preted as a molecu- 
lar gas disk from the gas kinematics ([ Scovillc et alJll997t 
IDownes fc Solomonlll998t ISakamoto et alJll999f >. In dust 
emission, on the other hand, the extended component is 
weakly detected o r not detected at a signifi cant signal 
level. At 1.3 mm. IDownes fc Solomon! (|1998[ ) suggested 
that the flux of 55 ± 11 mJy can be a ttributed to the 
flux fr om the extended component, but Sak amoto et al.l 
(1999) did not detect any signifi cant emission from the 
extended component. At 860 /xm. ISakamoto et all (|2008l ) 
attributed about 25% of the total flux density emission 
from outside the two nuclei. 

Assuming that all of the missing flux of our 435 /xm 
continuum data is from the extended component, we de- 
rived the dust properties of this extended component us- 
ing the x 2 fitting mentioned above using the data be- 
tween 1.3 mm and 435 fim. We adopt for the size of 
this component the same size that we adopted for the 
fitting of the total flux density, which is l'/94 x l'/28 
(|Scoville et al.l Il997h . The fitting result is shown in 
Fig. [2] with a dash-dotted line, and the derived values 
are T d ~ 38 K, (3 ~ 2.4, and the v c ~ 1500 GHz 
(~ 200 fim). These values are roughly consistent with 
the values derived by iGonzalez- Alfonso et al.l (|2004h us- 
ing ISO LWS data; they derived for the extended com- 
ponent T d - 50 K, v c ~ 3000 GHz (~ 100 urn), and a 
source size of ~ l'/8 — l'/9 (these values vary depending 
on their models) assuming (3 of 2.0. These results sug- 
gest that a significant amount (> 50%) of 435 /im flux 
is in the extended component as mentioned in Sect. 13.11 
contrary to the longer wavelengths where continuum is 
dominantly from the two nuclei. 

A recent SMA U/LIRG survey shows that many of the 
sample galaxies (9 out of 14) have large missing fluxes 
(typically 50 — 80%) in 880 /im continuum emission, even 
though many of the galaxies also show compact distribu- 
tions. This result suggests that many of ULIRGs have 
extended c ontinuum emission with very luminous com- 
pact cores (| Wilson et al.ll2008l ). similar to the results of 
Arp 220. 

4.3. Excitation Conditions of Molecular Gas in the Two 

Nuclei 
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TABLE 3 

Integrated 12 CO intensities and intensity ratios at each nucleus 



Source 


I( 12 CO 2-1) 
[K km s" 1 ] 


I( 12 CO 3-2) 
[K km s- 1 ] 


I( 12 CO 6-5) 
[K km s- 1 ] 


12 CO(6-5)/(2-l) 


12 CO(6-5)/(3-2) 


12 CO(3-2)/(2-l) 


East Nucleus 
West Nucleus 


5540 ± 830 
4460 ± 670 


5040 ± 760 
5140 ± 770 


2390 ± 720 
2840 ± 850 


0.43 ±0.14 
0.64 ±0.21 


0.47 ±0.16 
0.55 ±0.19 


0.91 ±0.19 
1.15 ±0.24 



Note. — The adopted spatial resolution is l'/S x 0''8 with P. A. of 129°. The 12 CO(6-5) data for the eastern and western 
nuc lei of Arp 220 ar e our data, and those for uv- matched 12 CO(2-l) and 12 CO(3-2) data arc from Sakamoto ct al 
and ISakamoto et al.l ll2008fl . 
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Fig. 9. — Rotational transition dependence of the 12 CO brightness temperatures (SEDs) of the double nucleus in Arp 220 and those 
of other galaxies. The horizontal axis is the upper rotational transition levels of the 12 CO lines and the vertical axis is the brightness 
temperatures of the 12 CO lines on an arbitrary scale. We fixed the Z ( 12 CO) / (dv / dr) of 5 X 10" 5 (km s _1 pc -1 ) -1 . (a) Arp 220 12 CO 
SED of each nucleus overplotted with the LVG calculation results. Dashed lines are for two temperatures at a density of 10 3 ' 6 cm -3 . 
Dotted lines are for two densities at a temp erature of 50 K. (b ) 12 CO SED of Arp 220 and other galaxies. The multi-J 12 CO data o f the 
Galactic Center < 2.5°) are taken from IFixsen et al.l l|1999T ). and those of M82 and Mrk 231 are compiled bv lWeifi et al.l fl2005aT l and 
Papad opoulos et al.l (12007T) . r espec tively. The 12 C O data of BR 1 202-0 725 are taken from th e following papers: J=l-0: Riech eFs et al.l 
H2006I) . J=2-l: ICarilli et al.l J2002l) . J=4-3, J=7-6: IQmont et all [I1996T ), and J=5-4: IQhta et al] JT996T) . The best fitted molecular gas 
temperature curves, made from LVG calculations for a density of 10 36 cm -3 , are also overplotted for reference. 



We made 12 CO SEDs of the two nuclei using our 
12 CO(6-5) data with the interferometric 12 CO(2-l) and 
12 CO(3-2) data (|Sakamoto et al.ll999tl2008T) to study the 
12 CO excitation conditions. We first matched the short- 
est uv distance for all three data sets, and imaged at 
the same synthesized beam size of l'/3 x 0'.'8 with PA. of 
129°, which is the same spatial resolution as our 12 CO (6- 
5) image. The integrated intensities and line ratios of 
these three lines for each nucleus are listed in Table [3J 
and the relative 12 CO intensities of various transitions 
are shown in Fig. [5] for each nucleus. 

The decrease of 12 CO intensities toward higher- J in 
the western nucleus is smaller than that in the eastern 
nucleus, and indeed the intensity ratios are higher for 
the western nucleus than those for the eastern nucleus. 
These differences are, however, within observational er- 
rors. Thus the difference in the excitation conditions of 
molecular gas between the two nuclei is not significant. 

To discuss the excitation conditions more quantita- 
tively, we estimated the excitation conditions of the 
molecular gas u sing the large-velo c ity-gr adient (LVG) 
approximation (|Goldreich fe Kwanl Il974h . The col- 
lision rates for C O o f < 100 K were taken from 
Flower fc Launavl (fl985h and of > 250 K were from 
McKee et al] ([1982 ). In these calculations, we assume 
that all the 12 CO emission comes from the same region 
(i.e., one-zone model), and assume the 12 CO relative 
abundance over velocity gradient, Z( 12 CO) /(dv/dr), of 



two kinds of curves, one (dashed lines) is to see the tem- 
perature dependence (we fixed nn 2 of 10 36 cm" 3 ), and 
another (dotted lines) is to see the number density depen- 
dence (we fixed Tk of 50 K). We only plotted the curves 
close to upper or lower limits to show the possible ranges 
of the excitation conditions and the goodness of the fit- 
ting. Under the above LVG assumptions, it is estimated 
that the eastern nucleus has a molecular gas temperature 



of - 30-250 K, or a density of ~ 10 



3.5±0.2 



cm 3 . For the 



5 x 10 5 (km s pc 



western nucleus, we could only derive the lower limits, 
which are > 40 K for temperature and > 10 3 ' 5 cm" 3 for 
density. As mentioned above, the estimated molecular 
gas conditions for the two nuclei overlap, but the west- 
ern nucleus tends to have higher temperature or density 
than the eastern nucleus. This tendency is similar to 
that derived from the dust SED fitting (Sect. 14. ip . In- 
deed, the derived molecular gas and dust temperatures 
for both nuclei match well. This suggests that both the 
molecular gas and dust reside at similar regions and in 
thermal equilibrium. 

We also calculated the dependence of our results on 
Z( 12 CO) /(dv/dr). If we decrease Z( 12 CO) /(dv/dr) by 
an order of magnitude (i.e., decrease the 12 CO relative 
abundance or increase the velocity gradient or both; we 
fixed 71h 2 of 10 4 cm" 3 or Tk of 100 K), the temperature 
or the density increases by about a factor of a few in both 
nuclei. If we increase Z( 12 CO) /(dv/dr) by an order of 
magnitude (we fixed jjh 2 of 10 3 4 cm -3 or Tk of 30 K), the 
In Fig. [Ha), we overplotted temperature or the density decreases by about a factor of 
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a few. Therefore a small (within an order of magnitude) 
change in Z ( 12 CO) /(dv/dr) does not change the result 
significantly. 

Our 12 CO SEDs are compared with those of other 
galaxies in Fig. [9jb). We plotted the multi-J 12 CO 
intensities of the Galactic Center (normal and quies- 
cent galaxy), M82 (nearby starburst galaxy), Mrk 231 
(evolved ULIRG with an AGN at the nucleus), and BR 
1202-0725 (radio-quiet and CO bright high-z quasar at 
z of 4.69) together with the 12 CO SEDs for the Arp 220 
nuclei. The brightness temperatures of the Galactic Cen- 
ter quickly decrease with the increase of rotational tran- 
sitions, but those of BR 1202-0725, M82, and Mrk 231 
stay almost constant even at high-J transitions. We over- 
plotted the best fit temperature curves on each source for 
comparison (as shown in Fig. [9^, higher temperature can 
be replaced with higher density). The Galactic Center 
can be modeled well with low temperature (or low den- 
sity) conditions, but other galaxies are explained with 
higher temperatures (or higher densities). The two nu- 
clei of Arp 220 is similar to these higher temperature 
(density) galaxies, and different from the Galactic Cen- 
ter. This suggests that the molecular gas excitation con- 
dition in the double nucleus of Arp 220 is similar to these 
galaxies. Note that the 12 CO SED up to J=6-5 or 7-6 
data with the current accuracy are not enough to distin- 
guish the excitation conditions between these high tem- 
perature (density) galaxies, including the two nuclei of 
Arp 220. Higher accuracy or higher- J observations are 
needed to differentiate the excitation conditions of these 
galaxies. 

It is known that Arp 220, M82, and BR 1202-0725 have 
active star formation inside. The excitation conditions 
in the two nuclei of Arp 220 and in M82 are averaged 
over the central a few hundred pc, and that in BR 1202- 
0725 is averaged over several kpc. M82 has a gradient 
in the physical conditions from the center to the outer 
region (|Petitpas fe WilsonH ^OOl , and the physical condi- 
tions derived above is more similar to those of the center, 
where the starburst region exist. BR 1202-0725 has two 
sources, north and sou th, and the southe rn source may 
consist of two sources ()Carilli et al.ll200"2l) . probably in- 
teracting with each other. The similar excitation condi- 
tions of molecular gas regardless of the observed regions 
suggests that the observed molecular gas is biased to- 
ward the gas closely related to the star forming regions, 
and the effect of star forming activities to the exciting 
conditions of surrounding molecular gas is similar. 

The 12 CO SED study can also be a useful tool to 
search for AGN(s), since the nearby Seyfert galaxies ex- 
hibit strong enhancement of hi gher- J 12 CO lines toward 
AGNs (e.g.. iMatsushita et all l200l iHsieh et al.ll2008l) . 
Fig. Mb) exhibits, however, that the 12 CO SED of the 
AGN hosting ULIRG Mrk 231 does not display any 
higher-J enhancement, and 12 CO SED comparison be- 
tween Mrk 231 and the two nuclei of Arp 220 shows 
no clear difference. In addition, the comparison be- 
tween Mrk 231, starburst dominated galaxies M82 and 
BR 1202-0725 also show no clear difference. We there- 
fore could not find any evidence of an AGN in Arp 220 
from this 12 CO SED study. These results suggest that 
the AGN contribution to the surrounding molecular gas 
(at least for Mrk 231) is much smaller than the nearby 
Seyferts, possibly due to the smoothing effect by a larger 



(linear scale) beam or to a larger opacity effect toward 
the AGN. 

4.4. Molecular Gas Abundance Anomaly in the Central 
Region of Arp 220? 

As mentioned in Sect. 13. 5[ we obtained a very 
low 13 CO(2-l)/C 18 0(2-l) line intensity ratio of about 
unity. Recent SMA observations toward nearby ac- 
tive star forming galaxies (NGC 253, NGC 1365, 
and NGC 3256) show 13 C0 (2- 1)/C 1 8 0(2-1) ratios of 
- 4 (Sa kamoto etal] [2006aObL l2007h . The 13 CO(l- 
0)/C O(l-0) lin e ratios in nearby galaxies are ~ 4 
(jSage et al.lll991l ). similar to the J=2-l ratios. If both 
the 13 CO and the C 18 lines are optically thin, the 
13 CO/C 18 line ratios for J=2-l and J=l-0 are ex- 
pected to have almost the same values. Some of inter- 
ferometric 13 CO and C 18 observations of nearby galax- 
ies show 13 C0/C 18 ratio of about 2 (jMeier fc Turner! 
120041: IChou et al.l2007f) . but no t unity as in Arp 220 (note 
that some regions observed bv I Meier fc Turner! show the 
13 CO/C 18 ratios of about unity, but the S/Ns are low). 
The intensity ratio of about unity in Arp 220 is therefore 
unusual compared with those in other galaxies. 

The intensity ratio may be closely related to the abun- 
dance ratio; the intensity ratio is expected to be sim- 
ilar to the abundance ratio, if both lines are optically 
thin. The abundance ratio between 13 CO and C ls O in 
our Galaxy is 5.5 for the Solar System and 12.5 for the 
Galactic Center, and that in external galaxies is 3 — 5 
dHenkel fc MauersbergeHl993h . assuming [ 13 CO]/[C 18 0] 
= [ 13 C]/[ 12 C] x [ 16 0]/[ 18 0]. Indeed, the abundance ra- 
tios and the aforementioned intensity ratios for external 
galaxies are similar. The intensity ratio of about unity 
is unusual also from the abundance viewpoint. 

Here we discuss possible reasons for this low inten- 
sity ratio using the LVG calculations. We assume both 
13 CO and C 18 molecules are located at the same re- 
gion (one-zone model). Note that since the brightness 
temperatures are different between these two lines and 
the 12 CO(2-l) line (see Fig. [7J, it is evident that these 
two lines and the 12 CO (2-1) line emanate from differ- 
ent regions. We also assume that the [ 13 CQ1/ [C 18 Q1 
relative abundance ratio of 4 (Wa ng et al J 12004). Un- 
der this relative abundance ratio, both lines have to be 
optically thick for the line ratio to be unity. We calcu- 
lated assuming Z ( 13 CO) / (dv / dr) of 1 x 10~ 5 , 1 x 10~ 6 , 
and 1 x 10~ 7 (km s" 1 pc" 1 )- 1 . Z ( 13 CO) / (dv / dr) of 
1 x 10 -6 (km s _1 pc -1 )" 1 can be explaine d as the Galac- 
tic ab undances of [ 13 CO]/[H 2 ] = 1 x 10~ 6 (jSolomon et all 
Il979f ) with the velocity gradient of 1 km s _1 pc -1 . Other 
parameters are the same as in Sect. 14.31 

The calculation results are shown in Fig. \W[ In the 
case of Z( 13 CO) /(dv/dr) of 1 x 10~ 6 (km s" 1 pc" 1 )" 1 , 
a high density of nn 2 > 1 x 10 4 cm" 3 is needed even 
for Tk of 10 K, and about an order higher density is 
needed for 100 K to realize the 13 CO(2-l)/C 18 0(2-l) ra- 
tio of 1.0 ± 0.3. This is because both the 13 CO and 
C 18 lines easily become optically thin at lower- J with 
the increase of temperature, since the population moves 
to higher-J. To compensate this, the density, and there- 
fore the column density per unit velocity, A^( 13 CO or 
C w 0)/dv = Z( l3 CO or C 18 'O) / (dv / dr) xn Ha , also has 
to be high for both lines to be optically thick. If we 
increase Z( 13 CO)/ (dv/dr) by an order of magnitude, 
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Fig. 10.— The LVG calculation results for the 13 CO (2- 
l)/C ls O(2-l) ratio as a function of H2 number density and 
kinetic temperature. Solid, dashed, and dotted lines are the 
13 CO(2-l)/C ls O(2-l) ratios with Z ( 13 CO) / \dv / dr) of 1 x 10 -6 , 
1 X 10 -7 , and 1 X 10 — 5 (km s — 1 pc — 1 ) — x , respectively, under 
the [ 13 CO]/[C 18 0] relative abundance ratio of 4. Dot-dashed 
lines are the 13 CO(2-l)/C ls O(2-l) ratios with the [ 13 CO]/[C ls O] 
relative abundance ratio of 2 under Z ( 13 CO) / (dv / dr) of 1 X 
1CT 6 (km s' 1 pc- 1 )" 1 . 

the density decreases by about a factor of several at 
a certain temperature. This is because the increase of 
Z ( 13 CO) /(dv/dr) makes the line easier to be optically 
thick. The response will be opposite if we decrease 
Z( 13 CO) /(dv/dr) by an order of magnitude. 

In the case of a lower [ 13 CO]/[C 18 0] relative abun- 
dance ratio of 2 (half the abundance ratio used above 
with increasing [C 18 0]), the required density for the ra- 
tio of 1.0 ± 0.3 is about an order of magnitude lower at a 
certain temperature (Fig. [T0|) . This is because the C 18 
abundance increased from the above condition, the opac- 
ity and therefore the intensity of the C 18 line become 
similar to that of the 13 CO line. 

As shown above, the important parameters for the 
low 13 CO(2-l)/C 18 0(2-l) ratio arc (1) the molecular gas 
density and hence the column density per unit velocity, 
and (2) the molecular abundance. First, we discuss the 
gas density. The molecular gas density needs to be high 
(> 10 4 cm -3 ) for the ratio to be around unity. The 
molecular gas in Arp 220 indeed contains high density 
gas, which is supported by the observations of high crit- 
ical density molecules, such as HCN, HCO + o r CS (e.g., 
ISolomon etafl I1990L Il992t iGreve et al.l S008). On the 
other hand, there are many galaxies with the detections 
of these high critical density molecular lines, but almost 
no report of a 13 CO(2-l)/C 18 0(2-l) - 1 so far. One 
possibility of the difference between Arp 220 and other 
galaxies may be due to the large fraction of dense molec- 
ular gas. Our 13 CO(2-l) and C 18 0(2-l) images exhibit 
compact distribution around the double nucleus, and the 
molecular gas in A rp 220 is dominated by dense gas (e.g., 
IGreve et a l. 2008). If the dense gas is concentrated to- 
ward the double nucleus, most of the molecular gas to- 
ward the double nucleus can be dominated by dense gas, 
and this makes the column density high enough to result 
in the 13 CO(2-l)/C 18 0(2-l) of unity. 



Second, we discuss the molecular abundance. To re- 
alize the 13 CO(2-l)/C 18 0(2-l) intensity ratio of about 
unity with changing the abundance, two possibilities can 
be considered; one is the underabundance of the 13 CO 
molecule, and another is the overabundance of C 18 0. De- 
ficiency of 13 CO abundance is often suggested in merg- 
ing galaxies based on their larger 12 CO(2-l)/ 13 CO(2-l) 
ratios than those in starburst or Sey fert galaxies (e.g., 
lAalto et all 119911 : ICasoli et al.l ll992). But as is men- 
tioned in Sect. [331 the observed 12 CO(2-l)/ 13 CO(2- 
1) ratio of Arp 220 is not as extreme as those in the 
other merging galaxies, and not significantly different 
from those in starburst or Seyfert galaxies. In ad- 
dition, the possible reason for the 13 CO deficiency is 
the s elective photodisso ciation of the 13 CO molecules 
(e.g.. ICasoli et al.lll992j ). In this case, however, C 18 
molecules will be more affected by the selective pho- 
todissociation (Ivan Dishoeck fc Blacklll988t ICasoli et all 
Il992| ). hence this cannot be the cause. We therefore think 
that the underabundance of 13 CO is possible, but less 
likely. 

The overabundance of the C 18 molecule may be pos- 
sible to achieve under the circumstance of Arp 220. Mas- 
sive stars synthesize a large amount of the primary ele- 
ment, 12 C, at helium burning phase, and it goes into in- 
terste llar medium via supernova explosions (jCasoli et alJ 
I1992D. The 18 Q enrichme nt occurs also in massive stars 
(Hcn kel fc Mauersbergen 119931 : iSage et al.lll991l) . either 
Wolf-Rayet stars or ty pe II supernova e xplosions by 
partial helium burning (jAmari et al.lfl995l ). Since Arp 
220 is very active in star formation (Sect. [1]), both the 
12 C and 18 O enrichment due to the above mechanism 
can be realized. This can lead to the enrichment of 
the C 18 molecule. This possibility still needs to be 
studied quantitatively. Note that a recent molecular 
gas abundance study toward a young (several Gyr old) 
galaxy at z = 0.89 show low [ 13 CO]/[C 18 0] of 1.9 ± 0.2 
(|Muller et all [20061. This result also suggests that the 
low intensity ratio is related to an abundance anomaly 
during the young star formation epoch. 

5. CONCLUSIONS 

We observed the central region of Arp 220 in the 
12 CO(6-5), 12 CO(2-l), 13 CO(2-l), and C 18 0(2-l) lines, 
and 435 and 1.3 mm continuum simultaneously with 
the SMA. The two nuclei are clearly resolved in the 
435 fim image, and kinematically resolved in the 12 CO (6- 
5) image. 

For the double nucleus, we concluded as follows: 

• The difference of the peak intensities in our 435 fim 
image between the two nuclei is smaller than at 
longer wavelengths. From the dust SED fitting, 
the dust in the two nuclei is estimated to be opti- 
cally thick at 435 /im. The emissivity indices are 
estimated to be ~ 2.0 for the eastern nucleus and 
~ 1.0 for the western nucleus. This suggests that 
the dust properties, such as dust size distributions 
or dust compositions, are different between the two 
nuclei. 

• The 12 CO SEDs are similar between the two nu- 
clei with the western nucleus having higher upper 
limits in the excitation conditions than those in 
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the eastern nucleus. The CO SEDs for both nu- 
clei and that of M82 or BR 1202-0725 are similar, 
characterized with small intensity decreases up to 
J = 6-5 ( 12 CO(6-5)/(2-l) ratio of about 0.5). This 
suggests that the molecular gas in the two nuclei 
of Arp 220 has the similar excitation conditions as 
that in M82 or BR 1202-0725, which have a density 
of > 10 3 ' 3 cm~ 3 or a temperature of > 30 K. 

• We could not find any evidence of an AGN in Arp 
220 with the 12 CO SED study. There is no clear dif- 
ference in the 12 CO SEDs between the AGN host- 
ing ULIRG Mrk 231 and the double nucleus of Arp 
220 (and therefore M82 and BR 1202-0725). This 
suggests that the AGN heating is not important 
for molecular gas excitation conditions in the large 
scale (a few hundred to a few kpc scale). 

For the global characteristics of the molecular gas and 
dust in Arp 220, we concluded as follows: 

• Based on the large amount of missing flux in our 
data and other previously published evidence in 
molecular gas and dust, we suggest the existence 
of an extended component in the dust emission 
with its dust properties T<j ~ 38 K, (3 ~ 2.4, and 
v c ~ 200 fim. A recent SMA U/LIRG survey sug- 
gests that m any of U/LIRGs see m to have extended 
components (IWilson et al.|[200^ ) , so that having an 
extended dust component might be common. 

• The 12 CO(2-l) line spectrum shows stronger line 
intensity at the lower velocities than the higher ve- 



locities, but the spectra of the higher-J lines show 
the opposite, indicating that the higher velocity gas 
has higher density, higher temperature, or both, 
than the lower velocity component. 

• The intensities of the 13 CO(2-l) and C 18 0(2-l) 
lines are similar. This suggests that the molecu- 
lar gas in Arp 220 is dense enough to be optically 
thick in both lines, or the abundance of either line 
deviates from the values in other nearby galaxies. 
To explain the ratio with the density effect, Arp 
220 should have molecular gas largely dominated 
by dense gas, more than in other nearby galax- 
ies. Underabundance of 13 CO is possible, but it 
is less likely. Overabundance of C 18 is also pos- 
sible, considering the 12 C and 18 enrichment by 
high mass stars. 
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